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Recent advances in understanding of magnetohydrodynamic (MHD) turbu- 
lence call for revisions in the picture of cosmic ray transport. In this paper we use 
recently obtained scaling laws for MHD modes to obtain the scattering frequency 

! for cosmic rays. Using quasilinear theory we calculate gyroresonance with MHD 
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We provide calculations of cosmic ray scattering for various phases of interstellar 

O ■ 
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c3 ! cosmic ray evolution. We obtain the scattering rate and show that fast modes 

j> \ provide the dominant contribution to cosmic ray scattering for the typical inter- 

stellar conditions in spite of the fact that fast modes are subjected to damping. 
We determine how the efficiency of the scattering depends on the characteristics 
of ionized media, e.g. plasma (3. We calculate the range of energies for which the 
streaming instability is suppressed by the ambient MHD turbulence. 

Subject headings: acceleration of particles-cosmic rays-ISM: magnetic fields- 
MHD-scattering-turbulence 



modes (Alfvenic, slow and fast) and transit-time damping (TTD) by fast modes. 
We provide calculations of cosmic ray scattering for various phases of interstellar 
medium with realistic interstellar turbulence driving that is consistent with the 
velocity dispersions observed in diffuse gas. We account for the turbulence cutoff 
arising from both collisional and collisionless damping. We obtain analytical ex- 
pressions for diffusion coefficients that enter Fokker-Planck equation describing 



1. Introduction 

Most astrophysical systems, e.g. accretion disks, stellar winds, the interstellar medium 
(ISM) and intercluster medium are turbulent with an embedded magnetic field that influ- 
ences almost all of their properties. High conductivity of the astrophysical fluids makes the 
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magnetic fields "frozen in", and influence fluid motions. The coupled motion of magnetic 
field and conducting fluid holds the key to many astrophysical processes. 

The propagation of cosmic rays (CRs) is affected by their interaction with magnetic 
field. This field is turbulent and therefore, the resonant interaction of cosmic rays with 
MHD turbulence has been discussed by many authors as the principal mechanism to scatter 
and isotropize cosmic rays (Schlickeiser 2002). Although cosmic ray diffusion can happen 
while cosmic rays follow wandering magnetic fields (Jokipii 1966), the acceleration of cosmic 
rays requires efficient scattering. For instance, scattering of cosmic rays back into the shock 
is a vital component of the first order Fermi acceleration (see Longair 1997). 

While most investigations are restricted to Alfven modes propagating along an external 
magnetic field (the so-called slab model of Alfvenic turbulence), obliquely propagating MHD 
modes have been included in Fisk et al. (1974) and later studies ( Pryadko & Petrosian 1999). 
A more complex models were obtained by combining the results of the Reduced MHD with 
parallel slab-like modes have been also considered (Bieber et al. 1988). Here we attempt to 
use models that are motivated by the recent studies of MHD turbulence (Goldreich& Sridhar 
1995, see Cho, Lazarian & Yan 2002 for a review and references therein). 

A recent study (Lerche & Schlickeiser 2001) found a strong dependence of scattering 
on turbulence anisotropy. Therefore the calculations of CR scattering must be done using 
a realistic MHD turbulence model. An important attempt in this direction was carried out 
in Chandran (2000). However, only incompressible motions were considered. On the con- 
trary, ISM is highly compressible. Compressible MHD turbulence has been studied recently 
(see review by Cho & Lazarian 2003a and references therein). Schlickeiser & Miller (1998) 
addressed the scattering by fast modes. But they did not consider the damping, which is 
essential for fast modes. In this paper we discuss in detail various damping processes which 
can affect the fast modes. To characterize the turbulence we use the statistics of Alfvenic 
modes obtained in Cho, Lazarian & Vishniac (2002, henceforth CLV02) and compressible 
modes obtained in Cho & Lazarian (2002, henceforth CL02, 2003b, henceforth CL03). 

As MHD turbulence is an extremely complex process, all theoretical constructions should 
be tested thoroughly. Until very recently matching of observations with theoretical construc- 
tions used to be the only way of testing. Indeed, it is very dangerous to do MHD turbulence 
testing without fully 3D MHD simulations with a distinct inertial range. Theoretical ad- 
vances related to the anisotropy of MHD turbulence and its scalings (see Shebalin et al. 
1983, Higdon 1984, Montgomery, Brawn & Matthaeus 1987, Shebalin & Montgomery 1988, 
Zank & Matthaeus 1992) were mostly done in relation with the observations of fluctua- 
tions in outer heliosphere and solar wind. Computers allowed an alternative important way 
of dealing with the problem. While still presenting a limited inertial range, they allow to 
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control the input parameters well making it easier to test theoretical ideas. For instance, 
one of us involved in the compressible turbulence research, failed so far to observe in the 
simulations the characteristic features, e.g. the production of a substantial amount of the 
quasi-parallel modes that are predicted by some of earlier MHD turbulence theories and 
were used to justify the slab geometry approximation for cosmic ray scattering. In this re- 
spect, the observational evidence in favor of such modes could be interpreted as the result 
of angle-dependent damping that we touch upon in this paper. 

While actual MHD turbulence is and will be in the near future the subject of intensive 
research and scientific disputes, in the present paper we adopt a simple model of turbulence 
that seems to be consistent with present day numerical MHD simulations. 

Yan & Lazarian (2002, henceforth YL02) used recent advances in understanding of 
MHD turbulence (CL02) to describe cosmic ray propagation in the galactic halo. In this 
paper we undertake a detailed study of cosmic ray scattering rates in different idealized 
phases of the interstellar medium. In §2, we describe our model of turbulence, introduce the 
statistics of Alfvenic and compressible turbulence in various conditions, including both high 
P and low f3 cases. In §3, we describe the resonant interactions between the MHD modes 
and CRs. The scattering by Alfven modes and fast modes is presented in §4. We apply our 
results to different phases of ISM, including Galactic halo, warm ionized medium (WIM), hot 
ionized medium (HIM), and also partially ionized medium. In §5, we describe the streaming 
instability in the presence of background turbulence. Discussion of our results is provided in 
§6 while the summary is given in §7. 



2. MHD cascade and its damping 

2.1. Adopted model of MHD turbulence cascade 

A shortcoming of many earlier studies was that the crucial element for cosmic ray 
scattering, namely, MHD turbulence briefly describe the model of turbulent cascade adopted 
and refer the reader to the original papers for a detailed explanation. 

Understanding of the properties of MHD turbulence is a challenging problem dealt with 
by many researchers. This work is not the appropriate place to cite all the important papers 
that shaped the field. For instance, a very incomplete list of such references in a review by 
Cho, Lazarian & Vishniac (2003) contains more than a hundred entries. 

It is well known that linear MHD perturbations can be decomposed into Alfvenic, slow 
and fast waves with well-defined dispersion relations (see Alfven & Falthammar 1963). Im- 
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portant questions arise. Can the MHD perturbations that characterize turbulence be sep- 
arated into distinct modes? Can the linear modes be used for this purpose? Is it possible 
to talk about, for instance, an Alfvenic cascade in a compressible medium? The answer to 
these non trivial questions is critical to studying cosmic ray scattering. 

The separation into Alfven and pseudo-Alfven modes, which are the incompressible 
limit of Alfven modes, is an essential element of the Goldreich-Sridhar (1995, henceforth 
GS95) model of turbulence. The arguments there were provided in favor of Alfven modes 
developing a cascade on their own, while the pseudo-Alfven modes would passively follow 
the cascade. The drain of Alfvenic mode energy to pseudo-Alfven modes was shown to be 
marginal along the cascade. 

This model and the legitimacy of the separation into modes were tested successfully 
with incompressible MHD numerical simulations (Cho & Vishniac 2000; Maron & Goldreich 
2001, CLV02). This motivated a further study of the issue within the compressible medium. 

The separation of MHD perturbations in compressible media into different modes was 
discussed further in Lithwick & Goldreich 2001 and CL02. There arguments were provided 
why one does not expect the modes to get entangled into an inseparable mess even in spite 
of the fact that MHD turbulence is a highly non-linear phenomenon. 

The actual decomposition of MHD turbulence into Alfven, slow and fast modes was a 
challenge that was addressed in CL02, CL03. There a particular realization of turbulence 
with mean magnetic field comparable to the fluctuating magnetic field was studied. This 
setting is probably typical of the Galactic magnetic field and would allow us to use a sta- 
tistical procedure of decomposition in the Fourier space, where the basis of the Alfven, slow 
and fast perturbations was defined. For some particular cases, e. g. for the case of low (3 
medium, the procedure was benchmarked successfully 1 and therefore argued to be reliable 
(see CL03 for more details). 

We also want to stress that the numerical studies in CL02 and CL03 were not blind 
empirical studies. They were aimed to test theoretical expectations. The very fact that for 
all the tested cases the theory was confirmed, makes us confident about the theory. 

Unlike hydrodynamic turbulence, Alfvenic one is anisotropic, with eddies elongated 
along the magnetic field (see Higdon 1984, Shebalin et al 1983). On the intuitive level it can 
be explained as the result of the following fact: it is easier to mix the magnetic field lines 



1 In this case, the velocity of slow modes is nearly parallel to the local magnetic field, and therefore 
the decomposition in the real space as compared to the Fourier case decomposition is possible. The two 
decompositions provided identical results (CL03) 
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perpendicular to the direction of the magnetic field rather than to bend them. However, one 
cannot do mixing in the perpendicular direction to very small scales without affecting the 
parallel scales. This is probably the major difference between the adopted model of Alfvenic 
perturbations and the Reduced MHD (see Bieber et al. 1994). In the GS95 model as well 
as in its generalizations for compressible medium mixing motions induce the reductions of 
the scales of the parallel perturbations. 

The corresponding scaling can be easily obtained. For instance, calculations in CLV02 
prove that motions perpendicular to magnetic field lines are essentially hydrodynamic. As 
the result, energy transfer rate due to those motions is constant Ek ~ ^\l T ki where is the 
energy eddy turnover time ~ where k± is the perpendicular component of the wave 

vector k. The mixing motions couple to the wave-like motions parallel to magnetic field giving 
a critical balance condition, i.e., k±Vk ~ fcyV^, where k\\ is the parallel component of the wave 
vector k, Va is the Alfven speed 2 . From these arguments, the scale dependent anisotropy 
fcy oc fc^ 3 and a Kolmogorov-like spectrum for the perpendicular motions Vk oc A; -1 / 3 can be 
obtained (see Lazarian & Vishniac 1999). 

It was conjectured in Lithwick & Goldreich (2001) that GS95 scaling should be approxi- 
mately true for Alfven and slow modes in moderately compressible plasma. For magnetically 
dominated, the so-called low f3 plasma, CL02 showed that the coupling of Alfvenic and com- 
pressible modes is weak and that the Alfvenic and slow modes follow the GS95 spectrum. 
This is consistent with the analysis of HI velocity statistics (Lazarian & Pogosyan 2000, 
Stanimirovic & Lazarian 2001) as well as with the electron density statistics (see Arm- 
strong, Rickett & Spangler 1995). Calculations in CL03 demonstrated that fast modes are 
marginally affected by Alfven modes and follow acoustic cascade in both high and low (5 
medium. In what follows, we consider both Alfven modes and compressible modes and use 
the description of those modes obtained in CL02, CL03 to study CR scattering by MHD 
turbulence. 

The distribution of energy between compressible and incompressible modes depends, in 
general, on the way turbulence is driven. Both the scale of driving and how the driving is 
performed matter. CL02 and CL03 studied generation of compressible perturbations using 
random incompressible driving with the energy coming from the largest scales of the system. 
While in the interstellar medium the energy can be injected at different scales, the bulk of 
the injected energy comes from the large scales (see Mac Low 2002) which explains the power 
spectra of turbulence observed (see discussion in Lazarian & Cho 2003). If the injection of 
energy happens at velocity higher that the Alfven velocity up to the scale at which V\ ~ Va 



2 note that the linear dispersion relation is used for Alfven modes. 
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the turbulence has hydrodynamic character. At smaller scales Alfvenic, fast and slow modes 
cascades develop. At the large scale the turbulence is isotropic in this case for all modes. 
The opposite case of injection velocity less than Va results in the Alfven and slow cascade 
being anisotropic even at the injection scale, provided that the driving is isotropic. 

How much energy is injected into fast modes if the driving takes place at large scales? If 
energy is initially injected in terms of Alfvenic perturbation, then the energy in fast ~ SVj? 
and Alfven ~ 5Vj modes are related by (CL02), 



where C$ is the sound speed. This relation testifies that at large scales incompressible driving 
can transfer an appreciable part of energy into fast modes. However, at smaller scales the 
drain of energy from Alfven to fast modes is marginal. Therefore the cascades evolve without 
much of cross talk. In high (3 medium, however, the velocity perturbation 5V at injection 
scale may be larger than Alfven speed. The injection speed of fast modes in this case can be 
approximated as SVf/SV ~ (5V/C S ) 2 - Naturally a more systematic study of different types 
of driving is required. In the absence of this, in what follows we assume that equal amounts 
of energy are transfered into fast and Alfven modes when driving is at large scales. 



In many earlier papers Alfvenic turbulence was considered by many authors as the 
default model of interstellar magnetic turbulence. This was partially motivated by the fact 
that unlike compressible modes, the Alfven ones are essentially free of damping in fully 
ionized medium (see Ginzburg 1961, Kulsrud & Pearce 1969) 3 . In the paper below we will 
show that compressible fast modes are particularly important for cosmic ray scattering. For 
them damping is essential. 

At small scales turbulence spectrum is altered by damping. Various processes can damp 
the MHD motions (see Appendix A for details). In partially ionized plasma, the ion-neutral 
collisions are the dominant damping process. In fully ionized plasma, there are basically two 
kinds of damping: collisional or collisionless damping. Their relative importance depends on 
the mean free path in the medium (Braginskii 1965), 



3 This picture contradicts to an erroneous assumption of strong coupling of compressible and incompress- 
ible MHD modes that still percolates the literature on turbulent star formation (see discussion in Lazarian 
& Cho 2003). However, little cross talk between the different astrophysical communities allowed these two 
different pictures to coexist peacefully. 





2.2. Plasma effects and damping of turbulence 
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Imfp = v th T = 6 x 10 n cm(T/8000K) 2 /(n/cm- 3 ). 



(2) 



If the wavelength is larger than the mean free path, viscous damping dominates. If, on 
the other hand, the wavelength is smaller than mean free path, then the plasma is in the 
collisionless regime and collisionless damping is dominant. 

To obtain the truncation scale, the damping time T^ 1 should be compared to the cas- 
cading time r k . As we mentioned earlier, the Alfvenic turbulence cascades over one eddy 
turn over time (k±v k )~ 1 ~ (/c||Va) _1 . The cascade of fast modes takes a bit longer: 



where V is the turbulence velocity at the injection scale, V p h is is the phase speed of fast 
modes and equal to Alfven and sound velocity for high and low [3 plasma, respectively 
(CL02). If the damping is faster than the turbulence cascade, the turbulence gets truncated. 
Otherwise, for the sake of simplicity, we ignore the damping and assume that the turbulence 
cascade is unaffected. As the transfer of energy between Alfven, slow and fast modes of 
MHD turbulence is suppressed at the scales less than the injection scale, we consider different 
components of MHD cascade independently. 

We get the cutoff scale k c by equating the damping rate and cascading rate r k T d ~ I. 
Then we check whether it is self-consistent by comparing the k c with the relevant scales, e.g., 
injection scale, mean free path and the ion gyro-scale. 

Damping is, in general, anisotropic, i.e., the linear damping (see Appendix A) depends 
on the angle 9 between the wave vector k and local direction of magnetic field B. Unless ran- 
domization of 9 is comparable to the cascading rate the damping scale gets angle-dependent. 
The angle 9 varies because of both the randomization of wave vector k and the wandering 
of magnetic field lines. Consider fast modes, that will be shown to be the most important 
for cosmic ray scattering. Their non-linear cascading can be characterized by interacting 
wave vectors k that are nearly collinear (see review by Cho, Lazarian & Vishniac 2003). The 
possible transversal deviation 8k can be estimated from the uncertainty condition 5ut cas ~ 1, 
where 5u ~ V p hSk(5k/k). Combining with Eq.(3), we can get 



r k = u/k 2 v 2 k 



(k/L)- 1 / 2 x V ph /V 2 , 



(3) 




(4) 



The field line wandering is mainly caused by shearing via Alfven modes. The deviation of 
field lines during one cascading period of fast modes can be estimated as: 
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B -\L) k ~ (fcL)V4 \V ) ■ (5j 

Thus the variation of the angle between k and B 59 ~ [l/kL) 1 ^ is marginal at large k. 
In the presence of anisotropic damping this results in anisotropic distribution of fast mode 
energy at large k. 

With this input at hand, it is possible to determine the turbulence damping scales for a 
given medium. For waves in fully ionized medium one should compare the wavelength and 
mean free path before determining which damping should be applied. However, when we 
deal with the turbulence cascade the situation is more complicated: If the mean free path 
is larger than the turbulence injection scale, we can simply apply the collisionless damping 
to the whole inertial range of turbulence. In general one should compare the viscous cut off 
and the mean free path. If the cutoff scale is larger than the mean free path, it shows that 
the turbulence is indeed cut off by the viscous damping. Otherwise, we neglect the viscous 
damping and just apply the collisionless damping to the turbulence below the mean free 
path 4 . By comparing different damping, we find the dominant damping processes for the 
idealized ISM phases (see Tablel). 



2.3. Statistics of fluctuations 

Within random-phase approximation, the correlation tensor in Fourier space is (see 
Schlickeiser & Achatz 1993) 



< ^(k)5*(k') > /B 2 = S(k - k')M^(k), 
< Vi (k)B*(k') > /V A B = 5(k - k')C^(k), 

< ^(k)^(k') > /Vl = 5(k - k')^-(k), (6) 

where Bi, Vi are respectively the magnetic and velocity perturbation associated with the tur- 
bulence. For the balanced cascade we consider, i.e., equal intensity of forward and backward 
modes, Cy(k) = 0. 

The isotropic tensor usually used in the literature is 

M l3 (k) = C {5 13 - hk 3 /k 2 }k- ll l\ (7) 



4 We shall show that for some angles between B and k the damping may result in turbulence transferring 
to collisionless regime only over a limit range of angles. 
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The normalization constant Cq here can be obtained if the energy input at the scale L 
is defined. Assuming equipartition w fc = J dk 3 Ylt=i MhBq/8tt ~ B 2 /8n, we get C = 
L~ 2 / 3 /127T. The normalization for the tensors below are obtained in the same way. 

The analytical fit to the anisotropic tensor for Alfven modes, obtained in CLV02 is, 



My(k) 
tfy(k) 



L-i/3 

67T 



-Iijk 



-10/3 



exp(-L 1/3 |A;|||/A; 



f) 



(8) 



where = {5^ — kikj/k 2 } is a 2D tensor in x — y plane which is perpendicular to the 
magnetic field, L is the injection scale, V is the velocity at the injection scale. Slow modes 
are passive and similar to Alfven modes. 

According to CL02, fast modes are isotropic and have one dimensional energy spectrum 
E{k) oc A;" 3 / 2 . In low (3 medium, the corresponding correlation is (YL02) 



£-1/2 

Hi i Jv 

8tt 3 



cos 2 9 
1 



(9) 



where 9 is the angle between k and B, ifj, = kikjjk\ is also a 2D tensor in x — y plane 5 . 
The factor cos 2 9 represents the projection as magnetic perturbations are perpendicular to 
k. 

In high j3 medium, fast modes in this regime are essentially sound modes compressing 
magnetic field (GS95, Lithwick & Goldreich 2001, CL03). The compression of magnetic field 
depends on plasma (5. The corresponding x-y components of the tensors are 



M tJ (k) 



£-1/2 

2tt 



sin 2 9H ij k- 7/2 



cos 2 9/(3 
1 



(10) 



The velocity perturbations in high f3 medium are radial, i.e., along k, thus we have the factor 
sin 2 # and also V\jC 2 = 2/(3 from the magnetic field being frozen cj<5B ~kx (v* x B). We 
use these statistics to calculate cosmic ray scattering arising from MHD turbulence. 



5 Here we give only the x,y components for the perturbation since only they are used in the following 
calculations. The complete 3D tensors are given in Appendix B. The tensor is not formally a tensor of an 
isotropic vector field because the vectors of fast waves are directed in a particular way in respect to vectors 
B and k. By definition, the velocity vectors of fast modes are in the B and k plane. This is why the tensors 
given here are different from eq.(7) 
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3. Interactions between turbulence and particles 

Basically there are two types of resonant interactions: gyroresonance acceleration and 
transit acceleration (henceforth TTD). The resonant condition is uo — knv/j, = nVt (n = 
0, ±1,2...), where u> is the wave frequency, Q = Qo/l is the gyrofrequency of relativistic 
particle, \x = cos£, where £ is the pitch angle of particles. TTD formally corresponds to 
n = and it requires compressible perturbations. 

We employ quasi-linear theory (QLT) to obtain our estimates. If mean magnetic field 
is larger than the fluctuations at the injection scale, we may say that the QLT treatment 
we employ defines parallel diffusion. Obtained by applying the QLT to the collisionless 
Boltzmann-Vlasov equation, the Fokker-Planck equation is generally used to describe the 
evolvement of the gyrophase-averaged distribution function /, 



dt 



d_ 

d/j, 



9f „ Of 

D — + D — 

w dfi w dp 



+ 



p2 Qp 



P 



r, df „ df 
D — + D — 

pp dp 



where p is the particle momentum. The Fokker-Planck coefficients D^, D^p, D pp are the fun- 
damental physical parameter for measuring the stochastic interactions, which are determined 
by the electromagnetic fluctuations (Schlickeiser & Achatzl993): 

Adopting the approach in (Schlickeiser & Achatz 1993) and taking into account only the 
dominant interaction at n — ±1, we can get the Fokker-Planck coefficients (see Appendix 
B), 



D 
D 



mi 



pp 




kj_vj_ 



i —)-4' k,v 



ft 



ft 



)) 



M nn (k) + M cc (k) 
\ i2<t> \ M^ £ (k) " 



+ e 



■i2<t> 



ft 



ft 



;n) 



where ( = 1 for Alfven modes and ( = k\\/k for fast modes, k m i n = L -1 , k maa 
corresponds to the dissipation scale, m = ^mn is the relativistic mass of the proton, v± is the 
particle's velocity component perpendicular to B , 4> = &rct&n(k y / k x ) , £,1Z = (x ±iy)/y/2 
represent left and right hand polarization. 
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The delta function 5(fc||i>|| — u> + nQ) approximation to real interaction is true when 
magnetic perturbations can be considered static 6 (Schlickeiser 1993). For cosmic rays, k\\V\\ ^> 
iv so that the resonant condition is essentially kuv/j,— nQ = 0. From this resonance condition, 
we know that the most important interaction occurs at k\\ = k\\^ res = Q/v\\. This is generally 
true except for small ji (or scattering near 90°). 



4. Scattering of cosmic rays 
4.1. Scattering by Alfvenic turbulence 



As we discussed in §2, Alfven modes are anisotropic, eddies are elongated along the 
magnetic field, i.e., k± > k\\. The scattering of CRs by Alfven modes is suppressed first 
because most turbulent energy goes to k± due to the anisotropy of the Alfvenic turbulence so 
that there is much less energy left in the resonance point k\\ tres = £l/v\\ ~ r^ 1 . Furthermore, 
k± fcy means k± ^> r^ 1 so that cosmic ray particles have to be interacting with lots of 
eddies in one gyro period. This random walk substantially decreases the scattering efficiency 
The scattering by Alfven modes was studied in YL02, (see Appendix B). In case that the 
pitch angle £ not close to 0, the analytical result is 



D 



ft' 



D 



pp 



f 2 - 5 /i 5 - 5 



r[6.5, kmaxk 



m 2 V] 



(12) 



where F[a, z] is the incomplete gamma function. The presence of this gamma function in 
our solution makes our results orders of magnitude larger than those 7 in Chandran (2000) 
for the most of energies considered. However, the scattering frequency, 



v = 2^/(1 -/i 2 ), (13) 

are still much smaller than the estimates for isotropic and slab model (see YL02). As the 
anisotropy of the Alfven modes is increasing with the decrease of scales, the interaction with 
Alfven modes becomes more efficient for higher energy cosmic rays. When the Larmor radius 
of the particle becomes comparable to the injection scale, which is likely to be true in the 



6 Cosmic rays have such high velocities that the slow variation of the magnetic field with time can be 
neglected. 

7 We compared our result with the resonant term as the nonresonant term is spurious as noted by Chandran 
(2000). 
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shock region as well as for very high energy cosmic rays in diffuse ISM, Alfven modes get 
important. 

It's worthwhile to mention the imbalanced cascade of Alfven modes (CLV02). Our 
basic assumption above was that Alfven modes were balanced, meaning that the energy 
of modes propagating one way was equal to that in opposite direction. In reality, many 
turbulence sources are localized so that the modes leaving the sources are more energetic 
than those coming toward the sources. The energy transfer in the imbalanced cascade 
occurs at a slower rate, and the Alfven modes behave more like waves. The scattering by 
these imbalanced Alfven modes could be more efficient. However, as the actual degree of 
anisotropy of imbalanced cascade is currently uncertain, and the process will be discussed 
elsewhere 8 . 



4.2. Scattering by fast modes 

The contribution from slow modes is no more than that by Alfven modes since the slow 
modes have the similar anisotropies and scalings. More promising are fast modes, which 
are isotropic (CL02). With fast modes there can be both gyroresonance and transit-time 
damping (TTD) (Schlickeiser & Miller 1998). 

TTD happens due to the resonant interaction with parallel magnetic mirror force. The 
advantage of TTD is that it doesn't have a distinct resonant scale associated with it. TTD 
is thus an alternative to scattering of low energy CRs which Larmor radii are below the 
damping scale of the fast modes. Moreover, we shall show later that TTD can contribute 
substantially to cosmic ray acceleration (also known as the second order Fermi acceleration). 
This can be crucial in some circumstances, e.g., for 7 ray burst (Lazarian et al. 2003), and 
acceleration of charged particles (Yan & Lazarian 2003). Different from gyroresonance, the 
resonance function of TTD is broadened even for CRs with small pitch angles. The formal 
resonance peak k\\/k = V p h/vn favors quasi-perpendicular modes. However, these quasi- 
perpendicular modes cannot form an effective mirror to confine CRs because the gradient 
of magnetic perturbations along the mean field direction VyB is small. Thus the resonance 
peak is weighted out and the Breit-Wigner-type resonance function should be adopted (see 
eq.B4 and Schlickeiser & Achatz 1993). 

Here we apply our analysis to the various phases of ISM. 



8 Prcliminary results by one of us show that the inertial range over which the degree of anisotropy of 
imbalanced turbulence is small is very limited. 



-13- 



TO A A 

ISM 
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C 5 (km/s) 
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91 


8.1 


7 


0.91 


0.35 


n(cm -3 ) 


10~ 3 


4 x 10- 3 


0.1 


0.4 


30 


200 


/ m / p (cm) 


4 x 10 19 


2 x 10 18 


6 x 10 12 


8 x 10 11 


3 x 10 6 


10 4 


L(pc) 


100 


100 


50 


50 


50 


50 


B(/xG) 


5 


2 


5 


5 


5 


15 




0.28 


3.5 


0.11 


0.33 


0.42 


0.046 


damping 


collisionless 


collisional 


collisional 


neutral-ion 


neutral-ion 


neutral-ion 



Table 1: The parameters of idealized ISM phases and relevant damping. The dominant 
damping mechanism for turbulence is given in the last line. HIM=hot ionized medium, 
CNM=cold neutral medium, WNM=warm neutral medium, WIM=warm ionized medium, 
DC=dark cloud. 
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Halo 

Recent observations (Beck 2001) suggest that the galactic halos are magnetically dom- 
inant, corresponding to low j3 medium. For low (5 medium, we use the tensors given by 
Eq.(9). For the halo, since damping at the scale of mean free path is slower than cascading, 
fast modes cascade further to collisionless regime and are subjected to collisionless damping. 
As we see from Eq.(A2), the damping increases with unless is close to n/2. This means 
that those fast modes in the two narrow cones with axis parallel and perpendicular to the 
magnetic field are less damped. For those quasi-parallel modes the argument for the Bessel 
function in Eq.(Bl) is k± tan£//c|| jres < 1 unless £ is close to 90°. So we can take advantage 
of the anisotropy of the damped fast modes and use the asymptotic form of Bessel function 
for small argument J n (x) ~ (x/2) n /n\ to obtain the corresponding analytical result for this 
case (see Appendix B): 
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The contribution from those quasi-perpendicular modes is 
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As we can see in Fig.l, the quasi-perpendicular modes are confined in a much narrower cone 

k kw 

than quasi-parallel modes, which indicates that -j^- 3> -^ ± - Thus AD UU ^ D ua . This 



'Uli- 



IS understandable because quasi-perpendicular modes have the similar anisotropy as Alfven 
modes, which has been shown negligible for CR scattering. 

From Eq.(14), we see the key factor is (k± >c / k\\ tres ) = tan# c . Using Eq.(A2) and the 
relation r k T L ~ 1, we see the critical angle C decreases with the energy of CRs tan# c > 
2{rrii/Tim e f3kL) l / /i oc r^ 4 . Therefore scattering frequency in the halo increases slowly with 
energy. Since the variation of during cascade 80 ~ 1/(A;L) 1//4 C at the corresponding 
scales, we assume the only effect of this variation is to reduce C by 80. Fig.2a. shows 
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the result corresponding to a pitch angle £ = 45° and we will adopt this value for the 
following calculations. The contribution from TTD can be obtained by evaluating Eq.(B4) 
numerically. The contribution is mostly from quasi-perpendicular propagating modes for 
which the collisionless damping is minimal (see Eq.A2) 9 . Limited by the variation of 9, the 
smallest scale the cascade can reach is a scale corresponding to k c (7r/2 — 59). Combining with 
Eq.(4) and (5), we can solve Eq.(A2) iteratively and get k c (n/2 — 59) ~ 4.3 x 10 _14 cm _1 . 
The contribution of TTD to scattering is smaller than that of gyroresonance (see Fig. 2a and 
3a) except for small /i. 

WIM 

According to the results in §2, the warm ionized medium (WIM) is also in low (3 but 
collisional regime. The mean free path l m f p = 6 x 10 12 cm (tablel), which corresponds to 
the resonant scale of CRs with energy ~ lOGeV. Thus For CRs harder than lOGeV the 
resonant modes with k\\ = Q/v\\ are subjected to collisional damping. Since the viscous 
damping increases with 9, we can apply Eq.(14) to the gyroresonance in WIM as well. By 
equating the damping rate Eq.(A4) with the cascading rate Eq.(3), we can obtain tan# c > 
Va/ (A: 3 / 4 / 1 / 4 ) -sy6p^/r7o) oc r^ 4 . Thus according to Eq.(14) the scattering frequency curve in 
WIM gets steeper (see Fig.2b). The angle variation 59 ~ l/(kL) 1 ' 4 < 9 C for CRs > lOGeV, 
so we adopt the same treatment as above, namely, subtracting 59 from 9 C . For lower energy 
CRs, the only available modes are those with k in a small cone as the residual of the collisional 
damping at large scales. Certainly, these modes are also affected by collisionless damping. 
However, the pitch angle of these modes is so small that the collisionless process is marginal. 
For scattering of CRs with low energy or large pitch angles, TTD is an alternative and can 
be evaluated by Eq.(B4) (see Fig. 2b & 3b). 

Hot Ionized Medium 

Hot ionized medium (HIM) is in high (3 regime and we adopt the tensors given by 
Eq.(lO). The mean free path of the HIM is of the order of parsec and thus the fast modes in 
this medium are subjected to collisional damping (see tablel). We use Eq.(A4) to acquire 
the formal damping scale as a function of 9. As shown in Fig.l), the peak around 9 ~ 
53° indicates that damping is the weakest here. However, since the angle variation 59 ~ 
1/(A;L) 1 / 4 ~ 20° is large at the scale of mean free path 2 x 10 18 cm, the cascade doesn't get 
though mean free path. Therefore, there is no gyroresonance interaction for moderate energy 
CRs. The only available interaction is TTD. The scattering frequency is ~ 1.2 x 10 _13 s _1 
for the energy range considered. 



9 Thus as a special case, 6 function gives a good approximation in halo except for small \i. 



-16- 



Partially ionized media 

Partially ionized media are similar to HIM in the sense that fast modes are severely 
damped by the ion-neutral damping. The cascade is cut off before the resonant scales of 
most of the CRs we consider. Indeed, for the parameters chosen, we find that gyroresonance 
only contribute to the transport of CRs of energy > 0.8TeV in DC and > ITeV for WNM and 
CNM. Therefore only TTD contributes to the scattering of moderate energy CRs in these 
media. As expected, the scattering is less efficient in partially ionized media. The scattering 
frequencies for moderate energy CRs are 2.5 x 10 _14 s _1 in WNM, ~ 4.0 x 10~ 15 s -1 in CNM 
and 5.7 x lO^V 1 in MC. 

From the above results, we see that the fast modes dominate CRs' scattering in spite 
of the damping. From Eq.(B4), we see the ratio of scattering and momentum diffusion 
rates p 2 D fJifl / D pp ~ l//i 2 , which means that the scattering and acceleration by TTD are 
comparable. However, for gyroresonance, we know that p 2 D^/ D pp ~ v 2 /(min{(3, 1}V^) 
(see Eq.(14)). Consequently for momentum diffusion (or acceleration) TTD is dominant . 

A special case is that the cosmic rays propagate nearly perpendicular to the magnetic 
field, so called the 90° scattering problem. It should be noted that with resonance broadening 
associated with the turbulence, the requirement v\\ = V p h/ cos 9 > V p h is relieved. The 
contribution from TTD thus becomes dominant for sufficiently small /i. However, since quasi- 
linear approximation is not accurate in this regime, proper calculation should be carried out 
with non-linear effects taken into account (Owen 1974, Goldstein 1976). 

It should be pointed out that all the results above are very much dependent on the 
plasma f3, which is somewhat uncertain. But the basic conclusion is definite: fast modes 
dominate cosmic ray scattering except in very high f3 medium where fast modes are severely 
damped. This demands a substantial revision of cosmic ray acceleration/propagation theo- 
ries, and many related problems should be revisited. 

5. Cosmic ray self confinement by streaming instability 

When cosmic rays stream at a velocity much larger than Alfven velocity, they can excite 
by gyroresonance MHD modes which in turn scatter cosmic rays back, thus increasing the 
amplitude of the resonant mode. This runaway process is known as streaming instability. 
It was claimed that the instability could provide confinement for cosmic rays with energy 
less than ~ 10 2 GeV (Cesarsky 1980). However, this was calculated in an ideal regime, 
namely, there was no background MHD turbulence. In other words, it was thought that the 
self-excited modes would not be appreciably damped in fully ionized gas. 
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Fig. 2. — Scattering frequency v given by Eq.(13) vs. the kinetic energy Ek of cosmic rays 
(a) in halo, (b) in WIM. The V lines refer to scattering by gyroresonance and solid lines are 
the corresponding analytical results given by Eq.(14). The dashed line are the results from 
Eq.(B4) for TTD. 




Fig. 3. — The scattering frequency v vs. cosine of pitch angle fi of lGeV CR (a) in halo, (b) in 
WIM. The 'x' lines refer to scattering by gyroresonance and solid lines are the corresponding 
analytical results given by Eq.(14). The dashed lines are the results from Eq.(B4) for TTD. 
In (a), 5 function and Breit-Weigner function match well except for the part marked by 
diamond symbols showing that 5 function underestimate the contribution (see text). 
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This is not true for turbulent medium, however. YL02 pointed out that the streaming 
instability is partially suppressed in the presence of background turbulence (see more in 
Lazarian, Cho & Yan 2003). More recently, detailed calculations of the streaming instability 
in the presence of background Alfvenic turbulence were presented in Farmer & Goldreich 
(2003), where the growth rate of the modes of wave number k was assumed to be (Longair 
1997) 



r(Jfe) = Q N{ - E \ -1 + ^=), (16a) 
n p V A 

where N(> E) is the number density of cosmic rays with energy > E which resonate with 
the wave, n p is the number density of charged particles in the medium. The number density 
of cosmic rays near the sun is JV(> E) ~ 2 x 10" 10 ( J E/GeV)~ L6 cm- 3 sr- 1 (Wentzel 1974). 

Consider interaction with fast modes, which were shown earlier to dominate scattering. 
Such an interaction happens at the rate ~ [k/ Vj^^Vph/V 2 (see Eq.(3)). By equating 
the growth rate Eq.((16a)) and the damping rate Eq.(3), we can find that the streaming 
instability is only applicable for particles with energy less than 

lmax ~ 1.5 x l^[n-\V py jV){Lvn Q /V 2 t"} 1 / 1 -\ (16b) 

which for HIM, provides ~ 20GeV if taking the injection speed to be V ~ 25km/s. 

One of the most vital cases for streaming instability is that of cosmic ray acceleration 
in strong MHD shocks. Such shocks produced by supernovae explosions scatter cosmic rays 
by the postshock turbulence and by preshock magnetic perturbations created by cosmic 
ray streaming. The perturbations of the magnetic field may be substantially larger than 
the regular magnetic field. In this situation a different expression for the growth rate is 
applicable (Ptuskin & Zirakashvili 2003): 

akell 3 

( } " cV A (Y max - (1 + a)^)(kr g0 ni + ( C) 

where e is the ratio of the pressure of CRs at the shock and the upstream momentum flux 
entering the shock front, U is the shock front speed, a — 4 is the spectrum index of CRs 
at the shock front, r g0 = c/Q , A tot is the dimensionless amplitude of the random field 
AL = 5B*/Bl 

Magnetic field itself is likely to be amplified through an inverse cascade of magnetic 
energy at which perturbations created at a particular k diffuse in k space to smaller k thus 
inducing inverse cascade. As the result the magnetic perturbations at smaller k get larger 
than the regular field. As the result, even if the instability is suppressed for the growth rate 
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given by eq. (16a) it gets efficient due to the increase of perturbations of magnetic field 
stemming from the inverse cascade. 

Whether or not the streaming instability is efficient in scattering accelerated cosmic rays 
back depends on whether the growth rate of the streaming instability is larger or smaller 
than the damping rate. The precise picture of the process depends on yet not completely 
clear details of the inverse cascade of magnetic field. If, however, we assume that the small 
scale driving provides at the scales of interest isotropic turbulence the nonlinear damping 
happens on the scale of one eddy turnover time. If assuming the shock front speed U is low 
enough, then the advection of the turbulence by shock flow is negligible and the perturbation 
Atot "C 1. In this case, we can get the maximum energy of particles accelerated in the shock 
by equating the growth rate of the instability Eq.(16c) to the damping rate due to turbulence 
cascade Eq.(3): 

MLeB )°*U* /(0 . 5+B) 

' max - y m o.5y2 c 2 > ■ y LX)a > 

From this we can estimate 7 maa; ~2x 10 7 (t/kyr)~ 9 / 4 in HIM. 

As shown in the Eq.(16a) and (16c), the growth rate depends on the CRs' density. In 
those regions where high energy particles are produced, e.g., shock fronts in ISM, 7 ray burst, 
SN, the streaming instability is more important. 



6. Discussion 

6.1. Applicability of our results 

We have attempted to apply a new scaling relations obtained for MHD turbulence to the 
problem of cosmic ray scattering taking into account the turbulence cutoff arising from both 
collisional and collisionless damping. We considered both gyroresonance and transit-time 
damping (TTD). 

In earlier papers it was frequently assumed that Alfvenic turbulence follows Kolmogorov 
spectrum while the problem of CR scattering was considered. This is incorrect because 
MHD turbulence is different from hydrodynamic turbulence if magnetic field is dynamically 
important. We appeal to the earlier research (see review by Cho & Lazarian 2003a) which 
showed that MHD turbulence can be decomposed into Alfven, slow and fast modes. They 
have different statistics (see CL02). The Alfven and slow modes follow GS95 scalings and 
show scale-dependent anisotropy. 

Fast modes, however, are promising as a means of scattering due to their isotropy. Even 
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in spite of the damping fast modes dominate scattering if turbulent energy is injected at large 
scales. We provided calculations for various phases of ISM, including galactic halo, WIM, 
HIM and partially ionized media. As the fast modes are subjected to different damping pro- 
cesses, the CR scattering varies. In galactic halo, the fast modes are subjected to collisionless 
damping. In WIM and HIM ion viscosity dominates the damping. All these damping in- 
crease with 6 and marginal for parallel modes. As the result, the fast modes at small scales 
become anisotropic with k\\ ^> k±. In this sense, fast modes have slab geometry but with 
less energies. These modes can efficiently scatter CRs through gyroresonance. TTD is more 
important for large pitch angle scattering and for the regime where fast modes are severely 
damped. We find that use of 5 function entails errors and resonance broadening ought to 
be taken into account. The formal resonance peak cos 8 = V p h/v\\ is smeared out because 
perturbations along the mean field VnB is small for the quasi-perpendicular modes. In all 
the cases, the scattering efficiency decreases with plasma (5 because the damping increases 
with (5. 

As every other model, the model adopted has its own limitations. For instance, within 
our model we did not consider wavelike component (see Montgomery & Turner 1981, Kin- 
ney & McWilliams 1998). Such parallel propagating fluctuations are observed in the solar 
wind (Matthaeus et al. 1990) and included in some of the existing models of cosmic ray 
transport (see Bieber et al. 1994). One could argue that such fluctuations could potentially 
be described by the weak turbulence model (see Galtier et al. 2000) or the weak turbulence 
model with imbalance (see a special case in Lithwick & Goldreich 2003). 

However, weak balanced turbulence is known to have a limited inertial range: turbulence 
get strong very soon (see discussions of this point in Cho, Lazarian & Vishniac 2003). A 
work by one of us also shows that the same is true for the weak imbalanced turbulence 10 . 
Therefore, while such type of turbulence could indeed be important near the sources of 
turbulence, its implications over vast spaces of interstellar gas are limited. 

The same is true about our model of energy injection. We assume that the turbulent 
energy is being injected at large scales. While observations of power spectra are consistent 
with such a picture, it is clear that some energy is being injected at small scales. If the 
intensity of fluctuations produced by this local injection is higher than those created by 
the global turbulent cascade, such regions can be important places for scattering of cosmic 
rays. If Alfvenic turbulence which provides rather weak scattering when injected from large 
scales (Chandran 2000, YL02) were the only mechanism for scattering of cosmic rays, such 



10 In addition the propagation of the imbalanced turbulence is limited in a compressible medium (Cho & 
Lazarian, in preparation) 
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localized regions of weak turbulence or intensive local injection of turbulent energy would 
be the dominant places of scattering. As we identified fast modes as an efficient source 
of scattering, in the zeroth approximation in our treatment we disregarded those localized 
sources of scattering. 

We disregarded shocks as sources of scattering. We believe that shocks are more impor- 
tant means of cosmic ray acceleration than scattering. To simplify our model further we also 
disregarded magnetic mirrors created by molecular clouds (Chandran 2001), which statistics 
is uncertain and seem to be important most if other scattering means fail. 

There are other simplifications that were adopted in the paper. For instance, our model 
disregarded complications associated with the large scale Alfven perturbations. For the range 
of cosmic rays that we are dealing the Alfven fluctuations are well within the inertial range 
and therefore various non-linearities are not important. Due to the same reasoning, we do 
not deal with the effects of "bendover" associated with the energy injection scale (see Zank 
et al 1998). 

Similarly, we do not deal with the gyroresonance scattering of the particles moving to 
magnetic field at angles close to 90 degrees. This allows us not to consider complex kinetic 
effects that are important for the small scale turbulence (compare Leamon et al. 1999). 
In addition, we did not consider effects of turbulence in partially ionized gas below the 
scale at which kinetic energy fluctuations are damped by the viscosity. This new regime of 
MHD turbulence arises in the partially ionized gas (Cho, Lazarian & Vishniac 2002) and is 
characterized by the high degree of magnetic field intermittency. While our simple estimates 
do not show that this regime is important for cosmic ray scattering, effects of magnetic 
intermittency on cosmic ray transport do require more work. 

How good is the decomposition when the mean magnetic field is not strong is the issue 
that must still be tested. The decomposition procedure used in CL02 and CL03 requires 
existence of the mean magnetic field. However, one may argue (see CL03) that when mean 
magnetic field is weak, locally the nonzero magnetic field of larger eddies acts as the mean 
magnetic field for the smaller ones. Therefore physically the results obtained for the par- 
ticular model could be argued to be valid even in the absence of the mean magnetic field. 
Testing this numerically is possible by taking volumes in the real space, defining there the 
local direction of the local mean field that changes from one volume to another (see Cho, 
Lazarian & Vishniac 2003 for a detailed discussion) and then using a decomposition in the 
Fourier space as this is done in CL02. For a present numerical studies this is extremely 
challenging, however. Nevertheless, the correspondence of properties of MHD turbulence 
with and without mean magnetic field reported in Cho & Lazarian (2003b) supports the 
theoretical expectations about the modes. One way or another, there are theoretical expec- 
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tations that have been tested for the case most important for scattering of Galactic cosmic 
rays. These were being used for the present study. 

Although scattering via fast modes is isotropic at large scales and becomes slab-type at 
small scales at which damping gets important, our model is radically different from earlier 
discussed 2D+slab and 3D+slab models (Bieber et al. 1994, Schlickeiser & Miller 1998). We 
used theoretically motivated models of turbulence that have been tested numerically. There- 
fore the transitions from one regime to another as well as the values of the intensity of 
perturbations have good justification within our approach. 

We adopted a particular set of parameters of ISM while doing the calculation. The basic 
conclusion, namely, fast modes dominate the cosmic ray scattering, will remain the same even 
if the parameters are altered. However, since the f3 value in ISM remains uncertain from 
observation, the damping scale k c is somewhat uncertain. And this determines the energy 
limit down to which gyroresonance is applicable. 

6.2. Comparison with recent attempts in the same direction 

An earlier attempt to study the effects of the GS95 scaling on the propagation of cosmic 
rays is made in Chandran (2000) and YL02 . These studies showed that Alfven modes are 
not be efficient for cosmic ray scattering. Moreover, YL02 identified fast modes as the driver 
for cosmic ray scattering. As fast modes are subjected to damping it was not clear whether 
fast mode induced scattering is always efficient. This problem was dealt with in the present 
paper, which provides detailed calculations of the expected scattering rate by fast modes for 
various phases of the ISM and takes into account both gyroresonance and TTD scattering. 

Another issue that this paper deals with is the suppression of the cosmic ray streaming 
instability by turbulence. This effect was first noticed in YL02 (see Lazarian, Cho & Yan 
(2003), Yan & Lazarian (2003) for more detail). The calculations of the effect were presented 
by Farmer & Goldreich (2004) for the weakly perturbed magnetic field. As the magnetic field 
may get strongly perturbed due to the inverse cascade of energy as the streaming instability 
develops first at smaller scales, in this paper we provide calculations of the suppression of 
the streaming instability for the strongly perturbed magnetic field. 

7. Summary 

In the paper above we characterized interaction of cosmic rays with balanced interstellar 
turbulence driven at a large scale. Our results can be summarized as follows: 
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1. Our calculations show that fast modes provide the dominant contribution to scat- 
tering of cosmic rays in different phases of interstellar medium provided that the turbulent 
energy is injected at large scales. This happens in spite of the fact that the fast modes are 
more subjected to damping compared to Alfven modes. 

2. As damping of fast modes depends on the angle between the magnetic field and the 
wave direction of propagation, we find that field wondering determined by Alfven modes 
affects the damping of fast modes. At small scales the anisotropy of fast mode damping 
makes gyroresonant scattering within slab approximation applicable. At larger scales where 
damping is negligible the isotropic gyroresonant scattering approximation is applicable. 

3. Transient time damping (TTD) provides an important means of how cosmic rays 
can be accelerated by fast modes. Use of 5 function resonance entails errors, and therefore, 
resonance broadening is essential for TTD. Our study shows that it is vital for low energy 
and large pitch angle scattering. And it dominates scattering of cosmic rays in HIM and the 
partially ionized interstellar gas where fast modes are severely damped. 

4. Streaming instability is subjected to non-linear damping due to the interaction of the 
emerging magnetic perturbations with the surrounding turbulence. The energy at which the 
streaming instability is suppressed depends on whether on the inverse cascade of magnetic 
energy as the instability gets more easily excited for low energy particles. 

We acknowledge valuable discussions with Jungyeon Cho, Steve Cowley, Don Cox, 
Randy Jokipii, Vladimir Mirnov, and Reinhard Schlickeiser. We thank the referee for helpful 
comments. This work is supported by the NSF grant AST0307869 ATM 0312282 and the 
NSF Center for Magnetic Self Organization in Laboratory and Astrophysical Plasmas. 



Below we summarize the damping processes that we consider in the paper. 
Ion-neutral damping 

In partially ionized medium, viscosity of neutrals provides damping (see LY02). If the 
mean free path for a neutral atom, l n , in a partially ionized gas with density n tot = n n + rii 
is much less than the size of the eddies under consideration, i.e. l n k <C 1, the damping time 



A. A. Damping of MHD turbulence 
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where v n is effective viscosity produced by neutrals 11 . The mean free path of a neutral atom 
l n is influenced both by collisions with neutrals and with ions. The rate at which neutrals 
collide with ions is proportional to the density of ions, while the rate at which neutrals 
collide with other neutrals is proportional to the density of neutrals. The drag coefficient for 
neutral-neutral collisions is ~ 1.7 x 10~ 10 T(K) a3 cm 3 s _1 (Spitzer 1978), while for neutral- 
ion collisions it is ~ (f r 0"m) ~ 1.9 x 10~ 9 cm 3 s _1 (Draine, Roberge & Dalgarno 1983). 
Thus collisions with other neutrals dominate for n;/n n less than ~ 0.09T 3 . Turbulent 
motions cascade down till the cascading time is of the order of td am p- The maximal damping 
corresponds to k^ 1 ~ l n . If the neutrals constitute less than approximately 5%, the cascade 
goes below l n and is damped at smaller scale (see below). 

Collisionless damping 

The nature of collisionless damping is closely related to the radiation of charged particles 
in magnetic field. Since the charged particles can emit plasma modes through acceleration 
(cyclotron radiation) and Cherenkov effect, they also absorb the radiation under the same 
condition (Ginzburg 1961). The thermal particles can be accelerated either by the paral- 
lel electric field which can also be called Landau damping or the magnetic mirror (TTD) 
associated with the comoving modes (or under the Cherenkov condition k\\v\\ = uj). The 
gyroresonance with thermal ions also causes the damping of those modes with frequency 
close to the ion-cyclotron frequency (Leamon et al. 1998), though it is irrelevant to the 
low frequency modes since we deal with GeV CRs here. The collisionless damping depends 
on the plasma (3 = P gas /Pmag and the propagation direction of the modes. In general, the 
damping increases with the plasma (3. And the damping is much more severe for fast modes 
than for Alfven modes. For instance, the Alfven modes are weakly damped even in a high 
(3 medium, where fast modes are strongly damped. The damping rate 7^ = r d x of the fast 
modes of frequency u for (3 <C 1 and 9 ~ 1 (Ginzburg 1961) is 

^ x/ttB sin 2 9 r / m e . m e . . 1 . . _ 

r ^ = ^-T^ rxL — exp ^_^) + 5exp(-- — , A2 

4 cos# ]/ m H m H (3cos 2 9 (3cos 2 9 

where m e is the electron mass. The exact expression for the damping of fast modes at small 
9 was obtained in Stepanov 12 (1958) 



11 The viscosity due to ion-ion collisions is typically small as ion motions are constrained by the magnetic 
field. 

12 We corrected a typo in the corresponding expression. 
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t l = ^-u9 2 x 1 + = J— exp — g — ). 

When /3 ^> 1, we obtain the damping rate as a function of 9 from Foote & Kulsrud 1979, 

sin 2 9 f 2cu 2 /a for fc < fli/pV A 
L ~ cos 3 X \ 2Q^//3 for fc > ^//JVk 

where fij is the ion gyrofrequency. 
Ion viscosity 

In a strong magnetic field (f^Tj ^> 1) the transport of transverse momentum is prohib- 
ited by the magnetic field. Thus transverse viscosity r]± is much smaller than longitudinal 
viscosity rjo = 0.96n(T/eV)Ti, ~ i] / (QiTi) 2 . The heat generated by this damping is 
Qvis — Vo(dv x /dx + dvy/dy — 2dv z /dz) 2 /3, where v x , v y , v z are the velocity components 
of the wave perturbation (Braginskii 1965). From the expression, we see that the viscous 
damping is not important unless there is compression. Therefore Alfven modes is marginally 
affected by the ion viscosity. 

While the damping due to compression along the magnetic fields can be easily under- 
stood, it is counterintuitive that the compression perpendicular to the magnetic also results 
in damping through longitudinal viscosity. Following Braginskii (1965), the damping of per- 
pendicular motion may be illustrated in the following way. To understand the physics of 
such a damping, consider fast modes in low f3 medium. In this case, the motions are pri- 
marily perpendicular to the magnetic field so that dv x /dx = n/n ~ B/B. The transverse 
energy of the ions increases because of the conservation of adiabatic invariant v\jB. If the 
rate of compression is faster than that of collisions, the ion distribution in the momentum 
space is bound to be distorted from the Maxwellian isotropic sphere to an oblate spheroid 
with the long axis perpendicular to the magnetic field. As a result, the transverse pres- 
sure gets greater than the longitudinal pressure by a factor r, ~ h/n, resulting in a stress 
~ pTifi/n ~ i]Qdv x /dx. The restoration of the equilibrium increases the entropy and causes 
the dissipation of energy with a damping rate T ion = k\r)o/6pi (Braginskii 1965). 

In high (3 medium, the velocity perturbations are radial as pointed out in §2. Thus 
Qvis = Vo(k x v x + kyVy — 2k z v z ) 2 /3 = r] k 2 v 2 (sm 2 9 — 2cos 2 6*) 2 /3. Dividing this by the 
total energy associated with the fast modes E k = PiV 2 , we can obtain the damping rate 
r io „ = A; 2 7 ? o(l-3cos 2 ^) 2 /6p i . 

All in all, 



(A3) 
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(k 2 ± rio/6pi for/3<l 
lon U^oll-Scos 2 ^) 2 ^ for/3»l 1 J 

Resistive damping 

In the paper we ignore the resistive damping because of the following reason. The 
resistivity is much smaller than the longitudinal ion viscosity for fast modes. For parallel 
propagating fast modes which are not subjected to the viscous damping from the longi- 
tudinal viscosity, the resistive damping scale can be obtained by equating the damping 
rate c 2 k 2 '/Sno ± with the cascading rate of the fast modes, where <tj_ ~ 10 13 (T/eV) 3 / 2 (see 
Kulsrud & Pearce 1969) is the conductivity perpendicular to the magnetic field: k~ e \ = 
1.8 x l0 7 ( y L/pc) 1 / 3 (T/10 4 K)~ 1 (VA/kms~ 1 )~ 2 / 3 cm. This scale turns out to be much smaller 
than the mean free path, where collisionless damping takes over. For Alfven modes, the 
resistive damping rate is T res = c 2 A;^/87rcr|| (the term associated with fey is negligible be- 
cause of anisotropy), where cry = 1.96cx_i_is the conductivity parallel to the magnetic field. 
By equating it with the cascading rate of Alfven modes, we can get the resistive damp- 
ing scale k~} s ~ L(c 2 /%na\\ VaL) 3 / 4 . Comparing with the proton Larmor radius, we find 
k-}Jr L = 3.6 x lQ- 3 (L/pc) 1 ^(T/WK)- 3 l 2 n l l 2 (V A /C s ) 1 ^. It's clear that in many astro- 
physical plasmas the resistive scale is less than proton Larmor radius. At this scale the 
energy can be transfered either to protons or to electron whistler modes. According to Cho 
& Lazarian (2004) whistler modes are even more anisotropic that the Alfven ones. Therefore 
their heating of protons is unlikely. In terms of cosmic ray scattering they are of marginal 
importance. 



B. B. Fokker-Planck coefficients 

In quasi-linear theory (QLT), the effect of MHD modes is studied by calculating the 
first order corrections to the particle orbit in the uniform magnetic field, and the ensemble- 
averaging over the statistical properties of the MHD modes (Jokipii 1966, Schlickeiser & 
Miller 1998). Obtained by applying the QLT to the collisionless Boltzmann-Vlasov equation, 
the Fokker-Planck equation is generally used to describe the evolvement of the gyrophase- 
averaged particle distribution. The Fokker-Planck coefficients D^, D^, D pp are the funda- 
mental physical parameter for measuring the stochastic interactions, which are determined 
by the electromagnetic fluctuations: 

< B a (k)B* p (k>) >= S(k - k')P a/3 (k) < B a (k)E* p (k') >= S(k - k')T a(S (k) 

< E a (k)B;(k') >= 5(k - k!)Q aP {k) < E a (k)E* (k') >= 5(k - k')R a p(k), 
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Adopting the approach in Schlickeiser & Achatz (1993), we can get the Fokker-Planck 
coefficients, 



D 
D 
D 



lip 



pp 



2B 2 



dk 3 



i ( i _i_ [±Va 
\ ±_l " < 



mc 



m A c z 



dte -i(k in -oj+nSl)t { j2 +i 



. k±V J_ . 



L RccQ*) J 



+ Jn+l{—Q—)Jn-l{- 



where |k„ 



A;,, 



= L , |k max | = k max corresponds to the dissipation scale, where 
( = 1 for Alfven modes and ( = k\\/k for fast modes, dk 3 = 2irk±dk±dk\\, m = 'ym H is 
the relativistic mass of the proton, v± is the particle's velocity component perpendicular to 
Bo, 4> = arctan(/c 2/ / k x ) , £,1Z = (x ± iy)/y/2 represent left and right hand polarization. For 
low frequency MHD modes, we have from Ohm's Law E(k) = — (l/c)v(k) x Bo. So we can 
express the electromagnetic fluctuations Tij,Rij in terms of correlation tensors Kij. For uo 
we use dispersion relations u = k»VA for Alfven modes and u = kV p h for fast modes. Those 
dispersion relations were used to decompose and study the evolutions of MHD modes in 
CL02. 



For gyroresonance, the dominant interaction is the resonance at n 
over r, we obtain 



±1. By integrating 
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where 6=1 for low (3 medium and i = sin 2 6 for high beta medium. 
For TTD, 
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Define the integrand 
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^iri in 
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r 2 + (fc||U|| - uj) 2 

e i2 *M nc (k 
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i 2 m 2 V 2 A 



(B3) 



M nn (k) + M cc (k) + e i2 *M nc (k) + e-^M cn (k) 
-CW(k) - CcM - e^C n M - e'^Ccn^) 
K nn {k) + K cc ik) + e i2 1>K nc (k) + e-^K cn (k) J 



/ 



G(k) = (J 2 (^) + J 2 (^))(^(k) + K cc {k)) - 2J 2 (^)J (^)(e^K nc (k) + e-^KMV 

T(k) = 2J 1 2 (^)((^(k) + K cc {k) + (e^K nc (k) + e- i2 *K cn (k)) (B4 

. We shall show below how they can be simplified in various cases to enable an analytical 
evaluation of the integral. The spherical components of the correlation tensors are obtained 
in the following. 

For Alfven modes, their tensors are proportional to 

sin 2 — cos sin 



cos sin cos 2 



For Alfven modes, kj_v±/Q ^> 1 because of the anisotropy. So we can use the asymptotics 
of the Bessel function J n (x) = A/27vnr cos(a; — mr/2 — 7r/4) when x — > oo. Thus from 
Eq.(8), we simplify the integrand in Eq.(ll), G(k) oc A;^ 10 ^ 3 exp(— L^-^lk^/k 2 / 3 ), the integral 
of which will be oc (| k\\ |LV 3 )- 6 - 5 r[6.5, k^J 3 \k\\ \L 1 '% Then from Eq.(ll) we can get the 
analytical result for the scattering by Alfven modes as given in Eq.(12). 

For fast modes, their magnetic correlations have such forms 

cos 2 9 cos 2 cos 2 9 cos sin cos 9 sin 9 cos 
cos 2 9 cos sin cos 2 9 sin 2 cos 9 sin 9 sin 
cos 9 sin 9 cos cos 9 sin sin sin 2 

Their velocity correlations are 

cos 2 cos sin 

Kij oc j cos sin sin 2 ] . (B5) 
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for (3 < 1 and K i:j oc kikj/k 2 for (3 > 1. 

In general, it's more difficult to solve the integral in Eq.(Bl) for the fast modes because 
they are isotropic. However, if taking into account damping, Bessel function can be evaluated 
using the zeroth order approximation. Thus from Eq.(ll,B,B5), we see for gyroresonance 
the integrand can be simplified as G(k) oc (J 2 - J ) 2 k- 7/2 oc k~ 7/2 . For TTD, T(k) oc 
J 2 (kj_vj_/Q)k~ 7 / 2 oc Ar 3 / 2 -y 2 sin 2 # oc £r 3 / 2 -y 2 (1 — V 2 h /v 2 ). Then we can solve the integral in 
Eq.(B4) analytically and get the scattering rate for fast modes. 
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